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ABSTRACT 

We present a new analysis of the properties of star-forming cores in the Perseus molecular 



(20051. Our 



cloud, identified in SCUBA 850 [im data originally presented by Hatchell et al. 
goal is to determine which core properties can be robustly identified and which 
the extraction technique. Four regions in the cloud are examined: NGC 1333, IC348/HH211, 
L1448 and L1455. We identify clumps of dust emission using two popular automated algo- 



ing the catalogues of Hatchell et al. ( 2007 1, we separate these clumps into starless. Class 
and Class I cores. Some trends are true for both populations: clumps become increas- 
ingly elongated over time; clumps are consistent with constant surface brightness objects (i.e. 
M °^ R^), with an average brightness « 4—10 times larger than the surrounding molecular 
cloud; the clump mass distribution (CMD) resembles the stellar intial mass function, with a 
slope a = -2.0±0.1 for clfind and a = -3.15 ±0.08 for gaussclumps, which sti-ad- 
dle the Salpeter value (a = —2.35). The mass at which the slope shallows (similar for both 
algorithms at M w 6M0) implies a star-forming efficiency of between 10 and 20 per cent. 
Other trends reported elsewhere depend critically on the clump-finding technique: we find 
protostellar clumps are both smaller (for GAUSSCLUMPS) and larger (for CLFIND) than their 
starless counterparts; the functional form, best-fitting to the CMD, is different for the two 
algorithms. The GAUSSCLUMPS CMD is best-fitted with a log-normal distribution, whereas 
a broken power law is best for CLFIND; the reported lack of massive starless cores in previ- 



ous studies (e.g. Hatch ell et al.|2007t [Hatchell & Fuller 2008) can be seen in the CLFIND but 



not the GAUSSCLUMPS data. Our approach, exploiting two extraction techniques, highlights 
similarities and differences between the clump populations, illustrating the caution that must 
be exercised when comparing results from different studies and interpreting the properties of 
samples of continuum cores. 

Key words: submillimetre - dust: extinction - stars: formation - stars: evolution - ISM: 
clouds - ISM: individual: Perseus. 



1 INTRODUCTION 

Stars form out of gas in the densest areas of molecular clouds. 
Early studies of molecular cloud structure (e.g. [Blitz & Shu|1980^ 
concluded that a division into discrete clumps of emission was the 
best description. Within such clumps, star formation occurs inside 
denser cores ^Myers & Benson|1983[ l. A growing consensus now 
maintains that molecular clouds have a scale-free structure gov- 
erned by turbulence (e.g. Elmegreen & Falgarone|[T996[ |Stutzki| 
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|et al.|[T998{ [Elmegreen & Scalo]|2004| l, with clumps only an ar- 
bitrary categorization. We may reconcile these viewpoints to some 
extent; the clump population has a power-law spectrum of mass im- 
plying the overall cloud has a similar distribution (e.g. [WilliamsTl 
[Blitz & McKee|[2000^ . Equally, the self-similarity of molecular 
clouds must break down where gravitational collapse becomes im- 
portant. 

This paper focusses on the decomposition of molecular clouds 
into clumps of emission. The utility of such a description de- 
pends on whether the located clumps accurately represent star- 
forming cores. Clumps are typically located in either dust contin- 
uum or spectral-line datasets, which both have their own advan- 



© 2009 RAS 



2 E. I. Curtis and J. S. Richer 



tages and drawbacks. Continuum observations select high-density 
cores, which tend to be self-gravitating, but associate more mass 
to clumps than exists in reality, since objects may be superposed 
along the line of sight (e.g. [Smith, Clark & Bonnell|2008, l. Clump- 
finding analyses on spectral-line cubes, with their added velocity 
dimension, should enable clumps to be more accurately assigned 
but are subject to larger uncertainties from the details of the radia- 
tive transfer than continuum emission. Ultimately, high signal-to- 
noise ratios are required in both cases for clumps to converge to- 
wards the underlying core population ( [Ballesteros-Paredes & Mac] 
|Low|2002| l. 

Nevertheless, many insights into star formation from clump 
populations remain compelling. |Motte, Andre & Neri| jl998^ 
first pointed out that the mass distribution of compact contin- 
uum clumps (the clump mass distribution; hereafter CMD) is re- 
markably similar to the initial mass function of stars (IMF). At 
high masses, the sl ope is consiste nt with a Salpeter power law 



(dN/dM oc M~^-^^ |Salpeter|l955| , considerably steeper than the 
T^^,„o,- ia,„ f^r nn r■^„,^r^c■ iriM lr\A/i A/f- 1 6 „ „ Blitz|l993 L At 



power law for CO clumps {AN /dM oc ' e.g. 
lower masses, the CMD slope becomes shallower, around the point 
that samples start to suffer from incompleteness. This suggests that 
the mass of a star is established at its earliest phases (see also |Alves^ 
[Lombardi & Lada|2007^ and would seem to rule out models where 
the shape of the IMF is formed later, through e.g. dynamical in- 
teractions jBate & Bonnell|[2005[ l. However, the mapping of the 
CMD on to the IMF is not straightforward and many evolutionary 
schemes would fit the present data (e.g. [Swift & Williams|2008| l. 

Many wide-field surveys are about to or have just begun, 
which will locate and characterize thousands of star-forming cores 
from continuum data, e.g. the James Clerk Maxwell Telescope 
(JCMT) and Herschel Gould Belt Surveys ( [Ward-Thompson et al.| 
[2007bl (Andre & Saraceno 2005 1. One of the key science drivers 
for these projects is the measurement of the CMD to high precision 
at low masses, into the brown-dwarf regime (< 0.08 M^). Thus, 
an examination of different source-finding techniques is particu- 
larly timely. Many different methods have been deployed in the 
literature to locate clumps in molecular-line and continuum data, 
from identifications by eye to automated algorithms. However, few 
studies compare the sources located with different techniques (a 
notable exception is 'Schneider & B rooks|2004) . In this paper, we 
closely compare the populations of continuum clumps found in the 
Perseus molecular cloud (hereafter simply Perseus) using the two 
most popular automated algorith ms, CLFIND ([W illiams, d e Geus &[ 
[BUtz|1994^ and GAUSSCLUMPS ( [Stutzki & GUst en 1990), to high- 
light their differences and biases. We re-examine SCUBA 850 \xm 
data, presented originally by [Hatchell et al.[ p005 ), in the four 
clusters of cores where we have complementary HARP spectral- 
line data jCurtis, Richer & Bucklel[2009) , namely: NGC 1333, 
IC348/HH21 1 (IC348 for short), L1448 and L1455. Although only 
investigating a sub-set of the SCUBA data limits our sample size, 
these are the sources whose kinematics we will investigate subse- 
quently and any CMDs from the spectral-line data will be directly 
comparable (Curtis & Richer, in prep.). The regions selected still 
encompass the majority of the SCUBA cores identified by |HatchelI| 
[et al . ( 2005 2007), 58 out of 103; so we can extrapolate any con- 
clusions to the entire cloud population. This work is divided into 
two principal parts. In Section|3] we identify clumps using the two 
algorithms before matching them to a catalogue of SCUBA cores, 
classified as protostellar or starless on the basis of their SEDs by 



[Hatchell et .^1 ( 120071 ) (hereafter [HOT} . Second, Sections |4] and |5] 
present an analysis of the physical properties of the cores before 
we summarize our conclusions in SectionjS] 



1.1 Nomenclature 

Many different terminologies have been used to describe the hier- 
archical structure in molecular clouds; we follow [Williams et al.[ 
(2000). Within molecular clouds, individual over-densities are 
termed dumps, these are the objects identified by automated al- 
gorithms and do not necessarily go on to form stars. Clumps may 
contain cores, which are the direct precursors of individual or mul- 
tiple stars. Every clump that does not contain an embedded object 
is referred to as starless. Of these starless cores, only a subset, the 
prestellar cores (formerly pre-protostellar cores, [Ward-Thompson[ 
[et al.|1994l ) will be gravitationally bound and thus go on to form 
stars. 



2 OBSERVATIONAL DATA 

We extracted fully calibrated and reduced SCUBA 850 |i.m maps 
across the four regions in Perseus we observed with HARP from 
the data presented by [Hatchell et al.[ j2005^ , where we refer the 
reader for details of the observations and processing. In short, the 
data were taken during 20 nights between 1999 and 2003. The sky 
opacity measured at 225 GHz varied from T225 = 0.039 to 0.080, 
with most of the data taken in good conditions (T225 ~ 0.05). 
The beam size is 14arcsec (0.017 pc at 250pc, our assumed dis- 
tance to Perseus) and the maps are sampled on a 3 arcsec grid. 
The rms noise level varies across the maps from typical values of 
CTrms = 26mJybeam"' in NGC 1333 and IC348 to higher values, 
32 and 46mJybeam~' in L1448 and L1455. Hereafter, to be con- 
sistent with the [HatcheIl et al. [ studies and therefore allowing a di- 
rect comparison of our source catalogues, we take the noise to be 
Crms = 35 mJy beam^', the typical value for the entire cloud. 



3 THE POPULATION OF DUST CLUMPS 

(Sub)millimetre continuum mapping selects young, cold cores (5- 
20 K) at the earliest stages of star formation. These largely self- 
gravitating cores are intermediate in properties between CO clumps 
and infrared young stellar objects (see [Ward-T hompson et al.[ 
|2007a, for a review). The dust's thermal emission is optically thin 
nearly everywhere at these wavelengths, accurately tracing density 
throughout the cloud. 



3.1 Algorithms 

We use CUPID I Berry et al. 2007 1^ to identify clumps of emission 
in the SCUBA 850 \im maps, exploiting the two most widely used 
algorithms: 



' Distributed as part of the Starlink software collection, see 
http: // starlink. jach.hawaii . edu. 
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(i) 'Clumpfind' or CLFIND ^Williams et al.||1994) , which con- 
tours the data at evenly-spaced levels, starting close to the peak 
value. All distinct contiguous areas are considered above each con- 
tour level in turn. If an area has no pixels previously assigned to 
a clump, it is marked as a new clump. Alternatively if there is a 
pre-existing clump assignment, that clump is extended to the next 
lowest contour level. When two or more different assignments exist 
in a area, the pixels are distributed using a friend-of-friends algo- 
rithm. 

(ii) GAUSSCLUMPS jStutzki & Gusten|ri990) , which uses a 
least-squares fitting procedure to break-down the data into a num- 
ber of Gaussian-shaped clumps, in a fashion similar to the CLEAN 
algorithms used in interferometric mapping. 

CLFIND is a partitioning method with each pixel only allocated to 
a single clump, whereas GAUSSCLUMPS permits any number of 
clumps to overlap at a single spatial position. This means GAUSS- 
CLUMPS is better at handling blended sources, although its output 
clumps necessarily conform to a strict Gaussian profile as opposed 
to the arbitrary shapes from CLFIND. 

The way each algorithm find clumps of emission is controlled 
by a number of tuning parameters. Our overall aim in this paper 
is to compare the populations of clumps found independently by 
CLFIND and GAUSSCLUMPS. We therefore try to find the same ob- 
jects with both algorithms, namely clumps larger than the beamsize 
with peaks greater than 4cTnns, but do not try to make the algorithms 
find exactly the same identifications clump-to-clump. In so doing, 
we follow the majority of previous studies and select parameters 
for each algorithm that are considered optimal for rejecting spu- 
rious detections. In an alternative approach, one could try to pro- 
duce clump catalogues which were as similar as possible. How- 
ever, the standard algorithm outputs we supply are quite different 
because of the very different way clumps are found with the two 
algorithms. Given that we can never definitively know the true, un- 
derlying clump population, then selecting one algorithm population 
to make more similar to the other would be completely arbitrary a 
priori. Indeed, it is only through detailed analysis, such as the one 
presented in this paper, or Monte-Carlo simulations that such a de- 
cision regarding algorithm selection could be made. 



3.1.1 CLFIND 

The most important parameters to tune CLFIND's output are the 
value of the lowest contour level {T\q.^) and the contour level spac- 
ing (AT). Ar of 2arms is considered optimal for rejecting clumps 
otherwise fitted to noise spikes ^Williams et al.|1994| l. T\g^ essen- 
tially sets the clump peak threshold (riow + ^T). After some ex- 
perimentation we settled on the optimum parameters to find plau- 
sible clumps: Ar and Tio^ both 2(7,n,s, so clumps must have peaks 
4CTi-nis. Additionally, the clumps had to be larger than the beam 
size and contain at least 7 pixels. We find 85 clumps in this man- 
ner which are depicted across the four regions in Fig. [T| The full 
catalogue of CLFIND detections is listed in Tab. |A1| 



Table 1. Number of clumps found for the best configuration of both CUPID 
algorithms alongside the number of sources in|H07| 



Algorithm\Region L1455 L1448 IC348 NGC 1333 Total 





H07 




5 


7 


17 


29 


58 


CLFIND 


7 


16 


23 


39 


85 


GAUSSCLUMPS 


8 


25 


24 


65 


122 



widths smaller than the beam size (14 arcsec) were rejected. In to- 
tal, 122 clumps were located in the four regions as shown in Fig.|2] 
Tab.|A2|provides the full the list of GAUSSCLUMPS identifications. 



3.1.3 Catalogue differences 

Although we identify clumps in the same data with the same algo- 
rithm, H07 find fewer clumps (the total number of clumps found in 
this work and H07^ is summarized in Tab.[TJ. H07 use the IDL imple- 
mentation of CLFIND on the SCUBA maps with an unsharp mask 
applied to remove structure on spatial scales > 2 arcmin. The most 
plausible explanation for this discrepancy is the different clump 
threshold, although there are other differences. In our catalogue, we 
include all CLFIND detections found with T\g^ = 2(7rms and Ar = 
2CTrtns- This ensures the detections all have peaks ^ 4(Trms- How- 
ever, HOT' only select their CLFIND detections (with Tjow = SfTmis 
and Ar = CJims) that have peaks ^ 5cTrms- Additionally, they select 
less bright sources, with peaks ^ ScTims that also have an identifi- 
cation by eye in their original paper ( [Hatchell et al.|2005] l or are 
Bolocam sources ( [Enoch et al.|2006| l. 



3.2 Distinguishing starless and protostellar clumps 

To classify the clumps as starless or protostellar, we looked for an 
associated source in the identifications of H07 The positional offset 
from each of our clumps' peak to the peak of its nearest [HOT] core 
in units of the observed clump diameter is plotted in Fig.[3](cf. sim- 
ilar analysis of |Enoch et al.|2008) . The diameter along each axis is 
twice the clump size, where the size is the standard deviation of the 
pixel coordinates about the centroid, weighted by the pixel values. 
We associate clumps with a core if they lie within one diameter. For 
CLFIND, 53 matching clumps (62 per cent of the total) were located 
within one diameter of the H07 source. All but two of these 53 are 
within half a diameter. Indeed most of our clumps lie directly on 
top of the |H07, sources so we can confidently classify them. The 
GAUSSCLUMPS results are similar with 52 clumps (43 per cent) 
matching to within one diameter. All but four of these lie within 
half a diameter. Three clumps are paired with two sources, as their 
emission is sufficiently large to encapsulate two |H07| cores. The 
separation of the |H07| cores and Gaussian clumps is significantly 
larger than the near perfect alignment of the CLFIND clumps. How- 
ever, the associations are still unambiguous and far from a random 
distribution. 



3.1.2 GAUSSCLUMPS 

For the GAUSSCLUMPS algorithm the most sensitive tuning param- 
eter was the threshold, the minimum value at which the routine 
would attempt a model fit (4cTrnis)- As for CLFIND, clumps with 



3.3 Dust mass estimates 

The mass of a clump, Mgjo, is proportional to its total 850 \im flux, 
5850, provided the emission is optically thin and its only source is 
thermal emission from the dust: 
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Figure 1. C lumps found with C LFIND in (a) NGC 1333, (b) IC348/HH211, (c) L1448 and (d) L1455. The regions selected are where we have complementary 
HARP data jCurtis et al.|2009} . The colour-scale is 850 yun SCUBA flux density from to 1000 mjy beam^' . The contours show the boundaries of the clumps 
and crosses mark the emission peaks, labelled with the clump number from the catalogue. 



My = 



d^Sy 
Bv{Td)Kv' 



(1) 



where Ky is the dust opacity, d the distance to the object and By the 
Planck function at Td, the dust temperature. Here the conversion 
becomes: 



^850 



0.70 5850 



exp 



X'850 
0.012cm2g-i 



V250pc 



Mr: 



(2) 



As discussed in [Curtis et al.| ( |2009) , we assume a constant 
distance of 250 pc to the entirety of the Perseus complex. Ini- 
tially, we adopt To = 12K everywhere. This temperature is con- 
sistent with the typical kinetic temperatures (T^ --^ 11 K) derived 



by [Rosolowsky et al.| l |2008) from ammonia observations towards 
193 dense core candidates in Perseus and dust temperature mea- 
surements at 5 arcmin resolution of the whole of the cloud I^Schnee] 
[eTaTIIIOOS) . [Schnee et al] (2008) find Perseus uniformly cool 
(Td ~12-20K) with a median = 15 K. When the sources are 
separated into classes by age, we assume a slightly higher tem- 
perature for those with an internal source of luminosity, i.e. pro- 
tostars (15 K), over starless cores (10 K). This is consistent with 
the average isothermal dust temperatures found in protostellar ra- 
diative transfer models (Shirley, Evans & R awlings''2002'; 'Young| 
et al.|[2003_l and c omplementary studies (|Hatchell & F uller 2008} 



Enoch et al.|2008| l. A final uncertainty in the mass comes from Ky 
We adopt a value, Ky = 0.012 cm^g^^, derived from the icy co- 
agulated grain model (n umber 5) of|Ossenkopf & Henning](|1994[l, 
which is consistent with |Hatchell et al.| ( |2005| ). |Kirk, Johnstone &] 
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Figure 2. Clumps found wi th GAUSSCLUMPS in (a) NGC 1333, (b) IC348/HH211, (c) L1448 and (d) L1455. The regions selected are where we have 
complementary HARP data jCurtis et al.|2009y The colour-scale is 850 \an SCUBA flux density from to lOOOmJybeam^'. The contours show the full- 
width half-maximum (FWHM) of the Gaussian model fit for each clump with crosses at the emission peak, labelled with the clump number from the catalogue. 



Di Francesco 1 2006 1, for example, prefer 0.02 cm^ ' , which re- 
suhs in our calculated masses being smaller by a factor of 1 .7 for a 
source of identical flux density. 



3.4 Completeness 

Completeness critically affects any CMD. At low masses, sensi- 
tivity limits cause a reduction in the number of clumps found and 
this could be confused with a turn-over in the mass spectrum. Fur- 
thermore, this completeness limit is a function of the clump mass, 
radius and profile. Massive, extended cores as well as all low-mass 
cores will be difficult to detect. Assuming constant surface bright- 
ness, a clump the size of the beam (radius 7 arcsec) must have at 
least a mass of ~0. IOM0 for detection. At the mean radius of the 
CLFIND clumps (21.4 arcsec, 0.026 pc) 0.95 is required, whilst 



an average sized (17.1 arcsec, 0.021 pc) GAUSSCLUIVIPS detection 
necessitates M 0.58 Mq. Uniform surface brightness is one of 
the simplest assumptions we could make for the clump profile. In 
reality, dust condensations are more centrally concentrated, with 
profiles well-fitted to Bonnor-Ebert spheres (e.g. [Johnstone et al.| 
pOOO; Kirk, Ward-Thompson & Andre||2005) . Higher concentra- 
tions ensure lower clump masses for a given radius and peak flux, 
so the uniform surface brightness estimate should be viewed as an 
upper limit. 



4 POPULATION CHARACTERISTICS 

In this section we explore the physical properties of the popula- 
tions, summarized in Tabs.|2]and|3] We examine how a core evolves 
over time by analysing the populations of different core ages, from 
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Figure 3. Distribution of the positional offsets (Aa,A5) between each 
clump peak found with CLFIND (top panels, (a) and (b)) or GAUSSCLUMPS 
(bottom panels, (c) and (d)) and the peak of its nearest'HOV'core by region. 
The distances are in units of the measured clump diameter The alignment 
is almost perfect for CLFIND and very good for GAUSSCLUMPS. 



prestellar to Class I sources, looking for significant differences in 
the distributions using Kolmogorov-Smirnov (K-S) tests. 

Overall, similar numbers of clumps are found for the two algo- 
rithms, apart from in L1448 and NGC 1333, where GAUSSCLUMPS 
finds considerably more. GAUSSCLUMPS fits a strict Gaussian func- 
tion, which is less flexible in a crowded, high-flux-density area 
with many clumps probably not following such a profile. Once one 
model has been subtracted there will still be ample flux in the resid- 
uals to continue to fit Gaussians, whereas a method such as CLFIND 
may distribute this flux to a smaller number of high-mass clumps 
(e.g. [Williams et al.|1994[|Smith et al.|2008) . 

4.1 Size 

In Figs. [4] and [3] we plot histograms of the clump deconvolved 
radius, ^dec^ formed as the geometric mean of the clump major 
and minor axis sizes (see Section [T2] l, each deconvolved with the 
beam half-width half-maximum (HWHM) size (^beam =7 arcsec or 
0.0085 pc): 



( ^dec.mai ■ 



maj ■^*^dec,min J 



Sizel^-i^beam •VSize2^--'^b 



1/2 



(3) 



The CLFIND population has a positively-skewed distribution with 
mean (i?dec) = (0.0244 ± 0.0014) pc |^ (see Tab. |2|. There are 
some minor differences region-to-region. IC348 and NGC 1333 
have similar distributions to the whole, the former skewed to 
slightly smaller sizes with (Rdec) = (0.020 ± 0.002) pc compared 
to NGC 1333's (0.0227 ±0.0016) pc. L1448 and L1455 have fewer 
clumps so their distributions are uncertain, but they do appear flatter 
with means of (0.027 ± 0.003) pc and (0.041 ± 0.006) pc respec- 
tively. K-S tests suggest the NGC 1333 population is significantly 
different from those in L1448 and L1455 with probabilities of being 
drawn from the same population of 7 and 2 per cent respectively. 
Similar results are found for IC348 compared to L1448 and L1455 
(5 and 2 per cent). 

The trends in the GAUSSCLUMPS data are similar (see Tab. 
[3](. The sizes are smaller than those for CLFIND by ~0.005 pc 
with an overall mean, (i^dec) = (0.0186 ± 0.0007) pc. The IC348 
and NGC 1333 clumps have similar means, (0.018 ± 0.002) and 
(0.0176±0.0009) pc respectively. The L1448 population is sfightly 
negatively skewed with the highest average, (Rdec) = (0.0225 ± 
0.0013) pc. K-S tests yield only 1, 3 and 4 per cent likelihoods that 
the L1448 population is drawn from the same one as NGC 1333, 
IC348 and L1455 respectively. Clumps in L1455 have changed 
greatly from their CLFIND distribution; they now have the small- 
est average, (i?dec) = (0.017 ±0.002) pc. 

For CLFIND the starless clumps are marginally smaller than 
the protostellar clumps, (Rdec) = (0.027 ± 0.003) compared to 
(0.0294 ± 0.0019) pc. The Class and Class I sources have nearly 
identical means ((0.0287 ± 0.0017) and (0.030 ± 0.003) pc) and 
similar distributions with the Class Os narrower. Intriguingly, with 
GAUSSCLUMPS the opposite trend is clear: on average starless 
cores are larger than Class sources, which in turn are larger than 



All the uncertainties listed throughout this paper are errors on the mean, 
formed from the sample standard deviations, (7. as a /\/N, where N is the 
number in the sample. 
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Table 2. Average clump properties for the CLFIND population. All the en'ors listed (cT) are errors on the mean not sample deviations. 



Population 






ARl" 


CTar 






M 






(pc) 


(pc) 






(10^3 cm--) 




(M.,.) 


(Ma) 


All 


0.0244 


0.0014 


1.50 


0.07 


1.40 


0.18 


6.4 


0.8 


NGC 1333 


0.0227 


0.0016 


1.44 


0.10 


1.7 


0.3 


6.8 


1.2 


IC348/HH211 


0.020 


0.002 


1.61 


0.17 


0.93 


0.18 


3.6 


0.7 


L1448 


0.027 


0.003 


1.50 


0.15 


1.5 


0.5 


9 


3 


L1455 


0.041 


0.006 


1.42 


0.10 


0.93 


0.16 


8 


2 


Starless 


0.027 


0.003 


1.46 


0.09 


1.00 


0.15 


6.0 


1.4 


Protostars 


0.0294 


0.0019 


1.42 


0.10 


2.6 


0.4 


10.6 


1.7 


Class 


0.0287 


0.0017 


1.37 


0.06 


3.4 


0.5 


12 


3 


Class I 


0.030 


0.003 


1.47 


0.19 


0.9 


0.2 


7.6 


1.8 



^ Deconvolved radius, s ee Se ction |4. 1 1 
* Axis ratio, see Section [4!2l 

Peak beam-averaged column density, see Section [431 
^ The entire Class and Class 1 protostellar population. 



Table 3. Average clump properties for the GAUSSCLUMPS population. All the en'ors listed (a) are errors on the mean not sample deviations. 



Population 






ARl" 


t^AR 


iV(H2)'^ 


On 


M 


t% 




(pc) 


(pc) 






(10^3 cm-2) 


(10^3 cm-2) 


(Me) 


(Me) 


All 


0.0186 


0.0007 


1.62 


0.07 


1.11 


0.12 


3.8 


0.3 


NGC 1333 


0.0176 


0.0009 


1.68 


0.11 


1.19 


0.18 


3.8 


0.5 


IC348/HH211 


0.018 


0.002 


1.69 


0.15 


0.85 


0.16 


3.0 


0.5 


L1448 


0.0225 


0.0013 


1.50 


0.10 


1.3 


0.3 


5.1 


0.8 


L1455 


0.017 


0.002 


1.25 


0.07 


0.74 


0.11 


2.3 


0.3 


Starless 


0.023 


0.002 


1.48 


0.10 


0.90 


0.09 


5.2 


0.7 


Protostars 


0.0182 


0.0011 


1.43 


0.07 


2.3 


0.3 


6.8 


0.9 


Class 


0.0190 


0.0012 


1.40 


0.10 


2.7 


0.3 


7.6 


1.1 


Class I 


0.017 


0.002 


1.49 


0.09 


0.72 


0.15 


5.4 


1.3 



^ Deconvolved radius, s ee Se ction |4.1| 
Axis ratio, see Section [4!2l 

Peak beam-averaged column density, see Section [431 
The entire Class and Class I protostellar population. 



Class I sources ((i?dec) = (0.023 ± 0.002), (0.0190 ±0.0012) and 
(0.017 ± 0.002) pc respectively). 

What do these results reveal about star-forming cores? In each 
region we might expect the clump sizes to vary according to the 
nature of the underlying population (i.e. protostellar or starless) 
and degree of clustering. Tab. |4] quantifies these effects, through 
the mean density of sources and the percentage of protostars. For 
CLHND the small clumps are in NGC 1333 and IC348 with larger 
clumps in L1455 and L1448. NGC 1333 is the most clustered envi- 
ronment so we would expect smaller clumps, since any measured 
size is limited by the distance to the nearest source. IC348 is an 
exception (low source density with small clumps), but its cores are 
significantly younger with many more designated starless that dom- 
inate the overall size statistics (CLFIND finds starless clumps are 
smaller than protostars). GAUSSCLUMPS should be better at sep- 
arating blended sources and indeed the clustering in the regions 
seems to have little bearing on the calculated sizes. The general 
trend in the population is for starless cores to be larger than proto- 
stars. This is the result reported by the Bolocam team ( [Enoch et al.| 
|2006[[2008) in their larger Perseus survey at 1.1 mm with a similar 



Gaussian fitting technique. Interestingly, they found the two types 
had similar sizes in Serpens and Ophiuchus. 

Clearly the type of clump analysis makes a big difference to 
the interpretation, a K-S test yields only a 1 per cent chance that 
the two distributions arise from the same population. In general 
CLFIND finds larger clumps than the tight spread of GAUSSCLUMPS 
radii. CLFIND produces arbitrary shaped clumps, rather than fitting 
a strict Gaussian profile. Therefore we might expect more emission 
at larger distances thus increasing the radii calculated. 

4.2 Shape 

The shapes of star-forming cores provide clues to their formation. 
Models with cores forming from turbulent flows predict random 
triaxial shapes with a slight preference for prolateness (jG ammiel 
et al. 2003] |Li et al.| [2004). If strong magnetic fields are present 
then cores are expected to be oblate ( ,Basu & Ciolek 2004, Ciolek] 
|& Basu|2006^ , which could also be caused by strong rotational mo- 
tion. Determining the exact 3-dimensional shape of a core from 
observations is impossible, instead statistical techniques have to 
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Table 4. Clump densities and protostellai' fractions. The protostellar fraction is the percentage of Class and we sources out of the total clumps associated 
with lH07l cores. 



Region 



Area 
/arcmin^ 



Number of clumps 



Mean clump density 
(arcmin^-) 



Percentage of protostars 
in identified sources 







CLFIND 


GAUSSCLUMPS 


CLFIND 


GAUSSCLUMPS 


CLFIND 


GAUSSCLUMPS 


NGC 1333 


190.0 


39 


65 


0.21 


0.34 


76.9 


72.4 


IC348 


198.6 


23 


24 


0.12 


0.12 


29.4 


31.3 


L1448 


111.1 


16 


25 


0.14 


0.23 


80.0 


71.4 


L1455 


127.1 


7 


8 


0.06 


0.06 


80.0 


80.0 




Figure 4. Distribution of clump radii deconvolved with the beam size, Rdco 
broken down by region for the CLFIND (top, (a) and (b)) and GAUSSCLUMPS 
(bottom, (c) and (d)) populations. 
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Figure 6. Distribution of clump axis ratios by broken down by region for the 
CLFIND (top, (a) and (b)) and GAUSSCLUMPS (bottom, (c) and (d)) popula- 
tions. The dot-dashed vertical line at a ratio of 1.2 separates clumps ,Enoch| 
|et al.|(2006} consider round from those that are elongated. 




0.01 0.02 O.Oa 0.04 



n.Ol 0-02 0-03 0-04 



Figure 5. Distribution of clump radii deconvolved with the beam size, Rdeo 
for the various core classifications of H07 for the CLFIND (top, (a) and (b)) 
and GAUSSCLUMPS (bottom, (c) and (d)) populations. 



be applied. Early axisymmetrical work indicated a preference for 
prolate cores (e.g. Myers et al. 1991 ) but later studies consistently 



favour triaxial, preferentially oblate shapes (e.g.|Jones, Basu & Du- 



binski 2001 ; Goodwi n, Ward-Thompson & Whitworth|2002t Tassis 
|2007||Tassis et al.|2009l i! 

A complete statistical analysis is beyond the scope of this 
paper, instead we examine the clump ellipticity in 2-dimensions, 
quantified via the axis ratio. The axis ratio is the ratio of the 
major to minor clump radii each deconvolved with the beam, 
^dec.maj/^decmin- We plot the distribution of axis ratios in Figs. 
6| and |7] whilst summarizing the data in Tabs. |2] and |3] |Enoch| 
et al.|(|2006 ) performed Monte Carlo simulations of a population of 



Gaussian sources to see the effects of noise and mapping technique 
on their Bolocam data in Perseus, which have a similar signal-to- 
noise ratio as ours. They conclude that a round underlying source 
can be distorted to reach a maximum axis ratio of 1.2. 

Region-to-region the distributions are similar, with K-S tests 
unable to reject the hypothesis of identical populations to any high 
level of significance. Most of the distributions have a positive skew 
with a large proportion of clumps having ratios close to unity. On 
average, the clumps are slightly elliptical with ratios of 1 .5 ± 0.07 
for CLFIND and 1.62±0.07 for GAUSSCLUMPS. Starless and pro- 
tostellar clumps have little difference in their degree of elonga- 
tion with ratios of 1.42 ±0.10 and 1.46 ±0.09 for CLFIND plus 
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Figure 7. Distribution of axis ratios by tlie source classifications of |H07| 
for tlie CLFIND (top, (a) and (b)) and GAUSSCLUMPS (bottom, (c) and (d)) 
populations. The dot-dashed vertical line at a ratio of 1.2 separates clumps 
[Enoch et al.H2006^ consider round from those that are elongated. 
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Figure 8. Distribution of clump peak beam-averaged column densities for 
the different regions and overall for the CLFIND (top, (a) and (b)) and 
GAUSSCLUMPS (bottom, (c) and (d)) populations. The densities were cal- 
culated assuming a constant dust temperature of 12 K. All regions exhibit 
similar behaviour apart from L1455 which has the fewest clumps. 



1.48 ±0.10 and 1.43 ±0.07 for GAUSSCLUMPS respectively. Class 
I sources appear marginally more elongated than Class 0, though 
the distribution is affected by a few extremely elliptical sources: 
average ratios of 1 .37 ± 0.06 and 1 .47 ± 0. 19 for Class and I ob- 
jects with CLFIND alongside 1.40 ±0.10 and 1.49 ±0.09 for the 
same with GAUSSCLUMPS. 

|ZerDovich| ( (T970l l and many others since have shown that a 
dynamically collapsing triaxial ellipsoid will collapse fastest along 
its shortest axis. Thus as collapse ensues axis ratios increase. This is 
more or less the picture we see here, with cores becoming slightly 
more elliptical as they evolve. Surprisingly, [Goodwin et ar]j2002^ 
found the opposite trend with starless cores more elliptical than 
protostellar ones in their more sophisticated statisti cal treatment of 
the deprojected shape of Jijina, Myers & Adams] ( |l999) 's N2H+ 
cores. This they tentatively propose as evidence of magnetic fields 
complicating simple gravitational collapse or of an outer protostel- 
lar envelope that is not collapsing at all. 



4.3 Column density 

The peak beam-averaged column density of hydrogen, A^(H2), de 
pends on the peak 850 jlva flux, fgjo, via: 



iV(H2 



(4) 



ilbeamM'«H2fv5v(7b) 

where flbeam is the solid angle subtended by the beam, mn^ the 
mass of a hydrogen molecule and /i = 1.4 the mean molecular 
weight per H2 molecule, assuming 5 H2 molecules for every He. 
For the values assumed in this paper: 



A'(H2 



: 7.62 X 10^2 

\iy beam 



exp I 1.4 
0.012cm2 g-i 



12K\ 



-1 



1 



-2 



(5) 



We take Tq to be 15 K for protostars, 10 K for starless agglom- 
erations and 12 K otherwise. Figs. [8] and |9] are histograms of the 



peak clump column density, whereas values for the individual pop- 
ulations are listed in Tabs. [2] and [3] These are necessarily popula- 
tion upper limits since lower density clumps are below our sensi- 
tivity threshold. For the first time the CLFIND and GAUSSCLUMPS 
distributions qualitatively resemble one another, with most clumps 
falling between 4 and 6.5 x lO^^cm^^, in the second bin. Their 
mean log column densities are also close, slightly higher than 
the most frequent bin at (1.40±0.18) x 10^^ cm"^ for CLFIND 
and (1.11±0.12) x lO^^cm"^ for GAUSSCLUMPS. The distribu- 
tions' regional variation does not seem to be significant, with only 
LI 455 (with the fewest clumps) showing markedly different be- 
haviour. The breakdown by age is more significant. Similar be- 
haviour is seen for the two algorithms, the protostars and starless 
cores have similar means, (2.6 ±0.4) x lO^-' and (1.00 ±0.15) x 
lO^-'cm"^ respectively for CLFIND alongside (2.3 ±0.3) x 10^^ 
and (0.90 ±0.09) x lO^-' cm"^ for GAUSSCLUMPS with the pro- 
tostellar sources having much wider distributions. When the pro- 
tostars are split into Class and I sources the lower densities are 
mainly Class I sources (means of (9±2) x 10^^ and (7.2 ± 1.5) x 
10^2 cmT^ for CLFIND and GAUSSCLUMPS) with the higher Class 
sources (means of (3.4 ±0.5) x lO^^ and (2.7 ±0.3) x 10^^ cm"^ 
for CLFIND and GAUSSCLUMPS). From the Class O/I definition, we 
have the simple expectation that the more evolved Class I cores 
have accreted more of their mass on to the central object. So there 
are lower densities in Class I envelopes, particularly given their 
smaller or comparable size to Class objects. 



4.4 Mass versus size relationship 

Fig.[TO]illustrates the mass-radius relation for the clumps. A least- 
squares fit to the whole population yields a straight line with a gra- 
dient of 1.7 ±0.1 for CLFIND and 1.2 ±0.2 for GAUSSCLUMPS. 
The slope is not well constrained; the constant surface bright- 
ness scaling, M oc R^, provides a good fit by eye. The scaling im- 
plies a surface brightness, E(H2) = 1920 MqPC^^ and £(H2) = 
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log(peak A'(H2)/cm-') 



log(i»akA-(Il2)/cm-») 




log(peaJ; A'(H2)/( 



Figure 9. Distribution of clump peak beam-averaged column densities for 
the clumps classified by age in the TIOT] catalogue for the CLFIND (top, (a) 
and (b)) and GAUSSCLUMPS (bottom, (c) and (d)) populations. We assume 
a dust temperature of 10 K and 15 K for the starless and protostellar clumps 
respectively. Protostars have a similar average density to the starless cores 
but a wider distribution. The Class sources are roughly the higher density 
half and Class Is the lower density half of the protostellar distribution. 



2470 M0 pc ^ for the CLFIND and GAUSSCLUMPS sources respec 
lively. These are around ten times higher than the av erage for 



molecular clouds as a whole, £(H2) = 17OM0pc"2 (|Solomon| 



|et al.|1987) and four tim es the largest Z found by I Solomon et al 
~500 I Heyer et al.|2009) . This implies there is a large density con- 
trast between the clumps and their surrounding medium. 

,Enoch et al.| l |2008} report distinct mass-radius relationships 
for protostellar and starless cores in Perseus, with the protostars 
deviating from the tight M correlation shown by their star- 
less counterparts. They calculate source sizes by fitting an ellipti- 
cal Gaussian (similar to GAUSSCLUMPS but at previously identi- 
fied peaks) and masses by integrating the flux within a fixed circu- 
lar aperture. Our results are slightly different depending on which 
algorithm is considered. The CLFIND relations are very close to 
M <x for whichever class of object, whereas for GAUSSCLUMPS, 
the scatter is greater around M '^Fp', with a somewhat greater range 
of masses for a narrower range of sizes in the protostars. This lat- 
ter observation loosely ties in with the simple suggestion of |Enoch| 
|et al.| |2008j that the present protostellar population in Perseus 
evolved from prestellar cores which decreased in size and density 
at constant mass until collapse and protostellar formation. We are 
only detecting the upper envelope of the mass-radius relationship, 
the bulk of the cores may exist below the completeness limit. How- 
ever, the small core numbers in-between the completeness limits 
for 10 and 15 K cores suggest that we are not missing very many 
protostellar objects. 

In summary, many of the previously reported properties of 
submillimetre dust clumps depend on the method of identification 
and care has to be taken to properly account for their biases. For in- 
stance, a simple evolutionary model with cores becoming smaller 
and denser with time would fit the GAUSSCLUMPS radii data with 
/?(Starless) > i?(Class 0) > i?(Class I). However, CLFIND data are 
contradictory; /^(Starless) < ./^(Class 0) ^ ./^(Class I). Other pat- 




Figure 10. Mass versus radius relation for the dust clumps in Perseus. The 
radius plotted, i?dec7 has been deconvolved with the beam (see Section [4~T) . 
The plots (all on the same scales) are for the CLFIND (left panels, (a), (c) 
and (e)) and GAUSSCLUMPS (right panels, (b), (d) and (f)) populations, dif- 
ferentiated by region (top panels, (a) and (b)) or by source type from the 
classifications of |H07| (the rest). Lines are drawn at various completeness 
limits for constant observed surface brightness i.e. with M <x R^, where R 
in this case refers to the observed clump radius (not deconvolved with the 
beam size). The limits correspond to 4oVnis detections above the noise for 
1 2 K clumps (dashed pink), 1 K starless cores (dashed black) and 1 5 K pro- 
tostars (dot-dashed green). Solid black lines mark the best-fitting straight 
lines of the form M c< R^ (i.e. constant surface brightness), either fitted to 
the whole population ((a) and (b)) or just the starless clumps (plots (c) to 
(f)). 



terns are more robust - the protostellar and starless populations are 
similarly elongated and the peak column densities of starless and 
protostellar clumps are close on average with the Class sources 
occupying the upper half of the protostars' distribution and the 
Class Is the lower. Larson's Law for constant surface brightness 
clumps, M 'X R^,m obeyed for both algorithms' clumps with more 
scatter in the protostellar population of GAUSSCLUMPS. 



5 THE CLUMP MASS DISTRIBUTION 

Every realistic model of star formation must match observations 
with the fundamental test being whether it reproduces the stellar 
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IMF. Interpretations of CMDs are complicated by a number of con- 
siderations, however appealing a direct one-to-one mapping may 
seem (see e.g. |Andre, Basu & Inutsuka|2009^ : 

(i) Form of the distribution. The differential form of the mass 
distribution, though easier to interpret, is inaccurate when the num- 
ber of sources is small (< 100), with a cumulative distribution pre- 
ferred ( |Reid & Wilson|2006| l. 

(ii) Dust properties. Core masses from continuum measure- 
ments are based on uncertain assumptions about the dust itself (i.e. 
opacity and temperature). 

(iii) Completeness. Sensitivity limits have the same effect as a 
real turn-over at lower masses, although in e.g. [Nutter & Ward-| 
|Thompson| ( |2007^ , a turn-over significantly higher than the com- 
pleteness limit is claimed. 

(iv) Nature of the objects. Many submillimetre cores without a 
central object will not form stars. Molecular-line observations can 
be useful in discriminating such unbound objects from their bound 
brethren (e.g. Curtis & Richer, in prep.). 

(v) Multiplicity. Some cores will form multiple objects, in- 
deed higher resolution observations often break cores down into 
smaller constituents. Furthermore, multiple systems can form after 
the prestellar phase during collapse (e.g.l Goodwin et al.|2007^ . 

(vi) Efficiency. Typically a constant star-forming efficiency, i.e. 
how effectively a core turns its mass into stellar material, is invoked 
to map the CMD on to the IMF. If this efficiency is mass depen- 
dent then the mapping can be more complicated ( [Swift & Williams] 
[2008| l. 

(vii) Timescales. [Clark, Klessen & Bonnell[ POOT} note that if 
core lifetimes depend on mass then the shape of the observed CMD 
will differ from the true, underlying distribution. 

We have already mentioned the striking similarity between the 
CMD and the stellar IMF and this section is a critical examination 
of this relationship. To begin, we examine the standard CMD for 
the SCUBA population, comparing it to previous work and looking 
at regional and age variations. Finally we will look at some of the 
distortions that prevent a simple mapping of the CMD on to the 
IMF. 



5.1 Form of the mass function 

For the initial analysis of the entire region, we mainly use a dif- 
ferential mass function, since its interpretation is more straightfor- 
ward. For clump numbers fewer than 100, where the arbitrariness 
of binning becomes important, the cumulative form is preferable 
( [Reid & Wilson|2006^ . We fit the differential CMD using the fol- 
lowing functional forms, the most simple being a single power law: 

^=AM« (6) 
dM 

with dA', the number of objects in a mass bin of width dM, a the 
power-law exponent and A a constant. Frequently, when a turn-over 
is observed (at break mass, Mb,eak)' a broken power law is fitted 
with two components: 



diV 
dM 



A(Mbreak)"'"*""'™^""" 



for M < Mbreak, 
for M ^ Mbreak 



(7) 



with (Zhigh and aiow the power law exponents above and below 
Mbreak respectively. Additionally, we fit a log-normal distribution 
to the differential CMD: 



AN 
dM 



Ma 



exp 



(lnM-/i)2 
2a2 



(8) 



where CT is the width, is related to the characteristic mass and A 
a normalization constant. 

It is normally assumed that the cumulative form of the CMD 
is also a power law differing in its exponent from a by unity: 

dA' 



N(>M) = / — dM = - - 
^ ./m dM a 



-M 



a+l 



1 



(9) 



for 0! < — 1 . In reality the sum is not to infinity but some fixed mass, 
Mmax, resulting from either an actual cut-off or sampling limit. Its 
shape is then modified: 

A^(> M) = rM"+' + -^^M„ 



a+l 



1 



1 



(10) 



This modification can introduce significant curvature at shallow 
slopes, which could easily be confused with the need to fit a broken 
power law. 



5.2 Differential mass function 

The differential mass distribution is presented in Fig. [TT| The 
clumps of both algorithms span a similar range of masses, with 
those from GAUSSCLUMPS extending to lower masses and cut- 
off at a smaller maximum mass: 0.3 to 47.8 Mg with an average, 
(M) = 6.4 M0 for CLFIND and 0.3 to 22.6 with (M) = 3.8 M© 
for GAUSSCLUMPS. Maximum-likelihood fits of a power law to 
the CLFIND data (above the sample average completeness limit 
as for all our fits), has an exponent, a — —1.5 ±0.1, and re- 
duced chi-squared, 2^ = 4.0. A broken power law provides a bet- 
ter fit ix^ = 0.9), with aiow = -0.6 ± 0.1 for M < 6.5 M0 and 
ajiigh = —2.0 ± 0. 1 for M > 6.5 M^. Additionally, a log-normal fit 
has parameters /i = 1.4 ±0.2 and CT = 1.2 ±0.1 with = 1.9. 
The GAUSSCLUMPS data are similarly poorly fitted by a single 
power law (a = —1.30 ±0.08, X^ = 3.7) but a broken power law 
is better (x^ = 0.5) , with aiow = -0.46 ±0.07 for M < 5.5 Mq 
and fthigh = -3. 15 ± 0.08 for M > 5.5 M©. The best-fitting distri- 
bution for the increased curvature is log-normal with parameters 
= 1 .0 ± 0. 1 and C7 = 0.97 ± 0.08, resulting in = 0.3. 

Our assumptions can affect the masses by up to a factor of 
four: e.g. dust temperature, opacity and the distance. Importantly, 
if these vary uniformly for every clump, the shape of the distribu- 
tion would be unchanged, merely translating it along the mass axis, 
altering the break mass. Non-uniform variations will certainly al- 
ter the shape. We will investigate temperature variations between 
clumps shortly but we completely ignore opacity and temperature 
changes across individual clumps which would require radiative 
transfer modelling. 

A CMD which closely resembles the IMF implies that the fi- 
nal masses of stars are set at the prestellar stage during the core 
formation process (e.g. [Padoan & Nordlund_2002) . In other mod- 
els of star formation e.g. competitive accretion, we would not ex- 
pect such a direct link with the earliest stages (e.g. [Bate & Bonneil] 
[2005[ l. A more realistic scenario might be intermediate between 
the two with dynamic interactions and turbulence playing a role 
( [Klessen "& Burkert.,2000^. The exac t shape of the local IMF is 
still rather uncertain (e.g. Scalo|2005| , with estimates of the power 
law slope from a = —2.3 to -2.8 at masses > 1 M©. For both of 
algorithms we find a break in the power law, significantly higher 
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Figure 11. Differential mass function of all the submillimetre clumps in 
Perseus across our four regions. The scales are the same for both plots 
as is the conversion to mass, with the standard assumptions including a 
constant clump temperature of 12 K. The error-bars represent Poisson un- 
certainties. Vertical lines (dashed grey) mark the completeness limits (0.95 
and 0.58 appropriate for CLFIND and GAUSSCLUMPS respectively) cal- 
culated for averaged-sized sources of each algorithm population (22 and 
IVarcsec in turn). Two further lines mark the IMF for single stars (thick 
dashed grey, |Kroupa|[200H and multiple systems (dotted grey, |Chabrier] 
[2005 1. Maximum likelihood fits to the mass function are plotted for vari- 
ous forms: a single power law (solid green), broken power law (solid red) 
and a log-normal distribution (dot-dashed blue). 



than the average sample completeness mass, although it should be 
stressed the sample is incomplete at every mass. At high masses 
we find a steeper power law, slightly lower than the Salpeter value 
for CLFIND and higher for GAUSSCLUMPS. A steeper slope im- 
plies a lack of high-mass objects (as is evident from the spread 
of masses). GAUSSCLUMPS probably misses clump mass at large 
distances from the clump centre, which can be found by CLFIND 
as it fits an arbitrary shape. The log-normal CMD fits are some- 
what wider than the jChabrier| ( |2005] ) IMF (1.2 and 0.97 compared 
to 0.55) and their characteristic masses are a factor of four or so 
larger (1.4 and 1.0 compared to 0.25 Mq). 

Assuming it is constant, we can estimate the efficiency of star 
formation, £, i.e. how much clump mass is converted into stel- 
lar, via the ratios of characteristic masses in the CMD and IMF. 
We compare the CMD break to the same in the IMF, (~0.5Mq, 
|Kroupa|20"0T] l yielding e = 0.08 and 0.09 for CLFIND and GAUSS- 
CLUMPS respectively. These are considerably smaller than the ac- 
cepted value of 0.3 for embedded clu sters ^Lada & Lada._2003j. 
Analysing all archived SCUBA data. Nutter & Ward-Thompson] 
l |2007t and fSimpson, Nutter & Ward-Thompson| ( |2008 | find smaller 
Mbreak of 2.4 Mq in Orion and 2.0 Mq in LI688 respectively. The 



different values suggest the efficiency depends on environment. 
The break mass can be heavily affected by any assumptions so we 
must exercise caution when making comparisons between differ- 
ent clouds. For example, the Perseus survey of [Enoch et al.|(2006| l 
finds the break at 2.5 M0, but if we scale our results to match their 
different assumptions about the dust opacity and clump tempera- 
tures, we find our break mass is reduced to 3.25 M© for CLFIND 
and similarly 2.25 M© for GAUSSCLUMPS. Correspondingly, e in- 
creases to 0.15 and 0.22 respectively and there is less difference 
between clouds. Nevertheless the lower values of £ than in clusters 
are probably because many cores will form multiple stars breaking 
the one-to-one mapping with the IMF (see [Goodwin et al.|2008[ 
[Swift & Wi lliams 2008 ). 

A number of further caveats complicate any simple interpreta- 
tion of the CMD. If higher mass clumps form more stars than lower 
mass ones, the CMD will be shallower than the IMF arising from 
it, as we found for CLFIND. Additionally, a clump's current mass 
may not be representative of the mass available for accretion over 
its lifetime or the star formation efficiency may vary with mass (see 
[Hatchell & Fuller|2008[ ). 

Other groups have measured the CMD in Perseus. The widest 
survey used Bolocam ( [Enoch et al.[[2006} which we plot along- 
side our results and the complementary SCUBA analysis of [Kirk[ 
[eTallpOOgl l in Fig. [T2I [Enoch et al.[ ( (2006 1 find 122 clumps, us- 
ing a bespoke algorithm with a more stringent detection limit, a 
5cTrms peak above the noise ((T,ms = 15mJybeam^') and larger 
beam (31 arcsec FWHM). Most sources in Perseus lie in five clus- 
ters, four of which we analysed, so our comparable numbers are not 
surprising. [Kirk et al.| ( [2006j l find considerably fewer sources, only 
58, even though they incorporated these data. They used their own 
data reduction technique, convolved to a degraded map resolution 
of 19.9 arcsec. The resulting maps had higher noise and they used 
the more stringent 5CTrms peak flux criterion with CLFIND, which 
probably explains the difference. 

[Enoch et al.| j2006[ l find a similar spread of masses (0.2 to 
25.6 M0), reaching out to lower masses due to their better sensitiv- 
ity. K^£rar] j2006j have a narrower range (0.2 to 9.8 M0), the 
absence of high mass clumps might arise as they have lost sensi- 
tivity to low density material with their higher noise and degraded 
resolution. Both groups adopted different assumptions; to match 
[Enoch et al.H2006| l and [Kirk et al.[ ( [2006| l, our data would need to 
be multiplied by 0.5 and 0.9 respectively. [Enoch et al.[ ( (2006[ l find 
their data is best fitted by a broken power law: = —1.3 ±0.3 
for M < 2.5 M0 and tthigh = -2.6 ± 0.3 for M > 2.5 M0. This is 
slightly steeper than for both of our algorithms at low masses and 
almost precisely intermediate between CLFIND and GAUSSCLUMPS 
at high. The jKirk et al.[ ( (2006| ) data are best fitted by a single power 
law that is quite shallow, a = — 1 .5 ± 0.2 but since there is a small 
number of clumps a cumulative mass distribution is more appropri- 
ate. 



5.3 Regional variations 

In Fig. [13] we plot cumulative CMDs for three of our regions in 
Perseus, L1455 has been excluded because it has a small number of 
clumps. We use the cumulative form as there are only a small num- 
ber of cores in each region. The outcomes of maximum-likelihood 
single and double power law fitting to the CMD with K-S tests to 
differentiate between the fits are presented in Tabs.[5]and|6] Clumps 
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Figure 12. Comparison of our differential mass functions, found with 
CLFIND (red) and GAUSSCLUMPS (blue) to previous work. We plot distri- 
butions derived by [Enoch et al.|(2006) (top panel, black) and |Kirk et al.| 
j2006) (bottom panel, green ), which conespond closely in shape to our 
mass functions. [Enoch et al.|j200"6) surveyed 7.5 deg^ of Perseus at 1 . 1 mm 
with Bolocam on the CSO, whereas [Kirk et aL|^200 6^ used the data pre- 
sented by [Hatchell et a l. 1 2005 1 combined with further independent SCUBA 
observations to yield a 3.5 deg^ survey. 



in L1448 are more massive than those in NGC 1333, which are 
more massive in turn than those in IC348. Returning to the vari- 
ous properties of the regions in Tab. |4] this implies that less clus- 
tered regions have more massive clumps (as we may expect from 
an analogous argument to that employed for the size) and regions 
that have higher protostellar fractions are more massive. The best- 
fitting models are typically single power laws with the maximum 
mass taken properly into account for CLFIND or broken power laws 
for GAUSSCLUMPS. The high-mass slope appears to be greater for 
IC348 than NGC 1333 or L1448, although the slope is very sensi- 
tive to the exact break mass and fitted model. This is not clear from 
the fits for GAUSSCLUMPS until the last points are looked at on the 
plot. Therefore there is a lack of high-mass objects in the IC348, 
which is likely to be caused by its young population, containing the 
largest proportion of starless cores. For instance, if starless cores 
are systematically colder than the assumed temperature (12 K), we 
would underestimate their masses. Additionally, their masses may 
be inherently different to protostars as we will explore shortly. The 
distribution in LI 448 varies markedly between the algorithms with 
the GAUSSCLUMPS slope being significantly steeper than the one 
for CLFIND. A shallow slope could be explained by its high pro- 
tostellar fraction, if they are over-estimated in mass by the 12 K 
assumption or are intrinsically heavier However, the steep GAUSS- 
CLUMPS slope for L1448 would bring us to the opposite conclu- 
sion. 
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Figure 13. Normalized cumulative mass distributions across three of our 
regions for the CLFIND (top, (a)) and GAUSSCLUMPS (bottom, (b)) popu- 
lations. The 4cT,-n]s completeness masses for a source of the average sample 
size are also plotted (dotted grey). 



Table 5. Best-fitting parameters for various models of the cumulative CMD 
split up by region. The models that are best-fitting to the data for each region 
are marked with an asterisk (evaluated using the K-S D values presented in 
Tah|6). 



Region 


1PL» 


IPLMM 


2PL 




\Modcl 


a 


a 


Kiow 


Chigh 


CLFIND: 










NGC 1333 


-L58±0.03 


-1.089 ±0.003* 


-1.47 ±0.02 


-2.39 ±0.04 


IC348 


-L79±0.05 


-1.15±0.02« 


-1.69±0.04 


-2.81 ±0.05 


LI 448 


-L46±0.06« 


no fit 


-1.32±0.01 


-2.39±0.06 


GAUSSCLUMPS 










NGC 1333 


-L68±0.02 


-1.43 ±0.04 


-1.60±0.02 


-2.55 ±0.04* 


1C348 


-1.56 ±0.08 


-l.O069±O.O0O3 


-1.41 ±0.04 


-2.43 ±0.10* 


LI448 


-1.72 ±0.07 


-1.23 ±0.02 


-I.13±0.01» 


-2.89 ±0.04 



" Single power law fit (IPL), as in Eq.[9](a everywhere corresponds to the 
gradient of an equivalent differential CMD). 

Single power law fit taking into account the finite upper limit of integration 
(i.e. maximum mass, IPLMM), as in Eq. [10[ 

Broken power law fit (2PL). 

5.4 Variations by source age 

An examination of how the CMD varies with source age can lead to 
crucial insights into the star formation process. In a trivial mapping 
on to the IMF, with a constant star-forming efficiency, the shape 
at the prestellar phase may distinguish between models where the 
mass is set by: (i) the initial fragmentation on to cores, (ii) compet- 
itive accretion or (iii) protostellar mechanisms i.e. outflows. 

In Fig. [14] we plot the cumulative CMD for protostars and 



© 2009 RAS, MNRAS 000,[7|fT8| 



14 E. I. Curtis and J. S. Richer 



Table 6. K-S D values for the various fits of the cumulative CMD by region. 
For the different model names see Tab.|5] An asterisk marks the best-fitting 
model for each region. 



Region & algorithm \Model 


IPL 


IPLMM 


2PL 


NGC 1333, CLFIND 


0.13 


0.06* 


0.17 


IC348, CLFIND 


0.09 


0.06* 


0.12 


L1448, CLFIND 


0.09* 




0.18 


NGC 1333, GAUSSCLUMPS 


0.13 


0.09 


0.08* 


IC348, GAUSSCLUMPS 


0.33 


0.12 


0.07* 


L1448, GAUSSCLUMPS 


0.26 


0.20* 


0.34 



Table 7. Best-fitting parameters for various models of the cumulative CMD 
split up by the source identifications from |H07] For the different model 
names see Tab.|5] The models that are best-fitting to the data for each source 
population are marked with an asterisk (evaluated using the K-S D values 
presented in Tab.jsJ. 



Region 


IPL 


IPLMM 


2PL 




\Model 


a 


a 


Clow 


"high 


CLFIND: 










Starless 


i.69±0.04 


-0.62±0.04» 


-1.57±0.04 


-5. 17 ±0.07 


Protostars - 


1.65 ±0.04 


-1.19±0.0I» 


-1.53±0.05 


-2.7I±0.05 


GAUSSCLUMPS 










Starless 


I.S2±0.04 


-0.68±0.05» 


-1.32±0.05 


-2.38 ±0.04 


Protostars 


l.J9±0.04 


-0.69±0.04» 


-1.48±0.05 


-3.71±0.04 



Starless cores with Fig.[T5]further breaking the protostars into Class 
and Is. The masses have been calculated in the usual way as- 
suming isothermal dust temperatures of 10 K for starless and 15 K 
for protostellar cores. As in Section [53] we fit models to the star- 
less and protostellar cumulative CMDs (see Tabs. [7] and [Sjl. For 
CLFIND, the protostellar CMD is wider, extending to higher and 
lower masses than the starless equivalent. The situation for GAtJSS- 
CLLFMPS is more complicated, with the protostellar distribution off- 
set to lower masses than the starless one. When the protostars are 
further divided, the Class sources have a flatter distribution, con- 
taining all the high mass objects, than the Class I or starless dis- 
tributions with CLFIND. The starless clumps are at lowest masses 
with the Class Is in the middle. The same holds for the GAtJSS- 
CLUMPS decomposition, except the Class I population extends to 
lowest masses and the starless population out to the highest. 

This is the part of the picture described in H07 and H atchell &| 
|FuIler| ( |2008p . For the same data analysed in a similar fashion with 
CLFIND, they find comparable results: Class sources are more 
massive than Class Is and there is an excess of Class protostars 
at high masses relative to starless cores - there is a lack of massive 



Table 8. K-S D values for the various fits of the cumulative CMD split up 
by source identifications from H07 For the different model names see Tab. 
[3] An asterisk mai'ks the best-fitting model for each category. 



Region & algorithm \Model 


IPL 


IPLMM 


2PL 


Starless, CLFIND 


0.14 


0.10* 


0.19 


Protostars, CLFIND 


0.35 


0.16* 


0.18 


Starless, GAUSSCLUMPS 


0.23 


0.06* 


0.09 


Protostars, GAUSSCLUMPS 


0.15 


0.06* 


0.11 
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Figure 14. Normalized cumulative mass distributions for protostars (green) 
and starless cores (black) as designated by |H07| for the CLFIND (top, (a)) and 
GAUSSCLUMPS (bottom, (b)) populations. The 4<7rins completeness masses 
for a source of the average sample size are also shown (dotted grey). The 
masses were calculated assuming a dust temperature of 10 K for the starless 
and 15 K for the protostellar sources. 



starless cores. It is reassuring that the Class sources are more mas- 
sive than Class Is. The original definition proposed sources to be: 
Class if the ratio of envelope to central object mass, Mgnv/Af* > 1 
and Class I if Menv/A^* < 1 i Andre, Ward-Thompson & Barsony 
[T993|. 



Our results are somewhat inconclusive. Single power law fits 
to the distributions derived by the two algorithms yield similar 
slopes for protostellar and starless sources. A power law fit, in- 
cluding the maximum mass term (Eq. |10^ , has a similar slope for 
both source types with GAtJSSCLUMPS and a shallower slope for 
the starless cores with CLFIND. This results because the maximum 
protostellar mass is considerably greater than the starless one, so 
the necessary curvature is already present and the slope can be shal- 
lower. A broken power law, however, has a steeper slope for starless 
over protostellar cores at the high mass end for CLFIND and the re- 
verse for GAUSSCLUMPS. 

All-in-all, the GAUSSCLUMPS distributions are very similar to 
one another. The Class I population is less massive than the Class 
and starless populations on average, which is again not surpris- 
ing given their definition. We expect more evolved sources to have 
accreted more mass and thereby have less massive dust envelopes. 
The CLFIND data support the view that the most massive objects 
are Class protostars and there is a paucity of massive starless 
cores. This could be explained by dust temperature variations or 
a mass-dependent evolution of the clumps, with massive prestel- 
lar clumps evolving rapidly into protostars, hence not as many are 
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Figure 15. Normalized cumulative mass distributions for starless cores 
(black). Class (blue) and Class I (purple) protostars as designated by |H07| 
for the CLnND (top, (a)) and G AUSSCLUMPS (bottom, (b)) populations. The 
4CTrms completeness masses for a source of the average sample size are also 
shown (dotted grey). The masses were calculated assuming a dust tempera- 
ture of 10 K for the starless and 15 K for the Class O/I sources. 
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Figure 16. Effect of varying clump temperatures on the mass distribu- 
tion. The original constant-temperature distributions (grey) have been over- 
laid with ones calculated using NH3 kinetic temperatures from the atlas of 
|Rosolowsky et al.|(2008) for the CLFIND (top, (a), red) and GAUSSCLUMPS 
(bottom, (b), blue) populations. The bars are standard Poisson errors on the 
bins and the 4f7n„s completeness masses for a source of average size are 
also shown (dashed grey). 



detected (see |Hatchell & Fuller|20d8t . However, if starless cores 
are not gravitationally bound, then they might not go on to form 
stars, so placing them in such an evolutionary sequence would be 
misleading. We examine their dynamic stability in a forthcoming 
paper (Curtis & Richer, in prep.). 



5.5 Effect of varying clump dust temperatures 

If the dust temperature varies between clumps then the shape of the 
mass distribution will change. Recently, Roso lowsky et af] j2008^ 
published NH3 kinetic temperatures, T^, towards 193 positions in 
Perseus using a beam size of 31arcsec. If the gas and dust are 
well-coupled, typically occurring at high densities, we would ex- 
pect Td ~ Tjf. Of these atlas positions, 50 coincide with our CLFIND 
identifications (59 per cent of the clumps, Tj^ = (12.9 ±0.3) K) and 
49 with GAUSSCLUMPS sources (40 per cent, = (12.6 ±0.3) K). 
In Fig. [16] we plot the mass distributions for the two algorithms 
using a varying dust temperature equal to where it was coinci- 
dent and 12 K where not. The distribution only changes marginally, 
indeed [Enoch et aLl ( |2008 ) found the same for a similar analysis 
on the Bolocam data - the deviation from a power law actually de- 
creased. This is not surprising since the average temperature of the 
clumps is very similar to our assumed dust temperature and their 
spread is small. 



5.6 Effect of clump lifetimes 

The observed mass distribution will not be representative of the true 
distribution if the lifetime of the cores is mass dependent, as we are 
more likely to detect long-lived cores than short-lived ones. |Andre| 
|et al.| ( |2007[ l suggest a weighted mass function, with each core 
weighted according to the inverse of it free-fall timescale ((%■) //ff, 
with (/ff) the average core free-fall time). Applying such weighting 
to the cores in Ophiuchus from |Motte et al.| ( |1998[ l, they find the 
break in the power law disappears and a single Salpeter-like power 
law is the best-fitting. In this section we attempt a similar weighting 
for our clumps. The free-fall timescale of a homogeneous sphere is: 

with p, the density in the clump, which we estimate from previ- 
ously derived quantities as shown. We find free-fall timescales ap- 
proximately in the range 10"* to 10^ yrs. The resultant weighted dis- 
tribution is shown in Fig.[T7] There is no apparent reduction in the 
flattening of the distribution at low masses. The reasons for this 
may be twofold. First, our estimates of the free-fall timescale may 
be in error. Other groups estimate densities using the mass and ra- 
dius inside a fixed aperture, whereas our radii are influenced by the 
whole of the clump. This will tend to reduce the average density 
of large radius clumps, i.e. large mass clumps, thereby increasing 
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Figure 17. Differential mass distribution weighted by the inverse of the 
clump free-fall timescale, (%)/%, with the original distributions in grey 
for the CLFIND (top, (a), red) and GAUSSCLUMPS (bottom, (b), blue) pop- 
ulations. The bars are standard Poisson errors on the bins and the 4t7niis 
completeness masses for a source of average size are also shown (dashed 
grey). 



their free-fall timescales and reducing their corresponding weight 
in the new distribution. This steepens the mass distribution at large 
masses and should therefore not affect the presence of a break. Sec- 
ond, there is, as other groups have noted, no relationship between 
mass and fff . If the break in the CMD is related to differential clump 
lifetimes then % should in some way be related to the mass. 



6 SUMMARY 

This paper detailed a new analysis of SCUBA 850 (im data towards 
NGC 1333, IC348/HH211, L1448 and L1455 in the Perseus molec- 
ular cloud. We used two clump-finding algorithms, CLFIND and 
GAUSSCLUMPS, to decompose the dust column density into clumps 
of emission. CLFIND located 85 clumps in total, whilst GAUSS- 
CLUMPS found 122. Once identified we divided the clumps into 
protostellar (Class and I) and starless sources by pairing them 
with classifications in the catalogues of |H07| 

This approach highlights trends in the clump properties that 
are highly or weakly dependent on the clump-finding algorithm. 
Conclusions which hold regardless of the technique include: 

• Shapes. Cores are slightly elongated with a tendency to be- 
come more elliptical over time as expected from collapse models. 

• Peak column densities. Protostars and starless cores have 
similar average peak column densities. The upper half of the pro- 
tostellar distribution arises from Class protostars and the lower 



from Class Is. This is expected from the class definitions and sup- 
ports the classifications of |H07| Unsurprisingly, the masses of Class 
clumps are also larger than Class Is. 

• Mass versus size. Clumps are consistent with the M re- 
lation, having average surface brightnesses some 4-10 times larger 
than typical giant molecular clouds. 

• CMD. The clump mass distribution resembles the IMF. The 
break mass is significantly higher than the average completeness 
mass and implies a star-forming efficiency of --^ 10-20 per cent, 
probably reduced relative to embedded clusters due to multiplicity. 

• Regional variations. The mass of clumps in each region 
follows: Meiu,np(L1448) > Me,u,np(NGC 1333) > M,iump(IC348). 
The high-mass slope of the CMD is steeper for IC348 than L1448 
with NGC 1333 intermediate. 

• CMD temperature and lifetime dependence. Using individ- 
ual temperatures for each clump or weighting according to their 
free-fall timescale does not change the shape of the CMD. 

Other inferences are dependent on the clump identification 
procedure: 

• Sizes. Protostars are either larger or smaller than prestellar 
cores with CLFIND or GAUSSCLUMPS respectively. 

• Functional form of the CMD. A log-normal CMD is best- 
fitting to the GAUSSCLUMPS data whereas a broken power law is 
best for CLFIND. In any case a broken-power-law fit has differ- 
ent slopes for each algorithm, with the GAUSSCLUMPS high-mass 
slope steeper at a = -3.15±0.08 than a = -2.0±0.I for CLFIND. 

• A lack of massive starless cores? A shortage of high-mass 
starless clumps relative to protostars was definitely present in the 
CLFIND data (see also |H07[ [Hatchell & F uller 20081. The same 
was not clear in the GAUSSCLUMPS data, where the overall spread 
in masses is smaller and the starless cores have larger sizes. 

GAUSSCLUMPS can only fit a strict shape so will not allocate 
flux to a clump at large distances from the peak unlike CLFIND. 
Therefore the corresponding GAUSSCLUMPS sizes are smaller. The 
clumps sizes in clustered regions such as NGC 1333 are smaller, 
as expected, for CLFIND but not any different for GAUSSCLUMPS, 
which is better at dealing with blended sources. The clump mass 
data are even more affected by selection effects. The slope of the 
starless and protostellar CMDs are very different for the two al- 
gorithms. Furthermore, the lack of high-mass starless cores, pre- 
viously observed by H07 and implied by the shallow prestellar 
mass function found by |Enoch et al.|p008^ is not apparent in the 
GAUSSCLUMPS data. This would imply that either the shortage is 
not real or GAUSSCLUMPS does not allocate enough flux to high- 
mass clumps. Indeed, protostars which have more irregular, dis- 
rupted envelopes due to e.g. outflows are perhaps worse affected 
by the regular profiles used by GAUSSCLUMPS. This would cause 
a lack of high-mass protostellar clumps, which might explain the 
GAUSSCLUMPS results. 

Fundamentally, we can never know the real population of 
clumps underlying any dataset. However, certain evolutionary 
trends must be present in any clump decomposition for it to be a be- 
lievable representation of this true population. One clear example 
of this is the evolution of core size over time. Every simple model 
of collapse naturally has cores decreasing in size with time. This 
trend is not apparent for the CLFIND clumps. In addition, particular 
CLFIND clumps incorporate material that seems equally likely to go 
on to form part of their final mass or that of their neighbours. There- 
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fore, treating the properties of a CLFIND clump as representative of 
gravitationally-bound material in a single core is probably highly 
misleading. Of course, we may need to tweak the CLFIND parame- 
ters, increasing the lowest contour level and decreasing the contour 
separation, but any method which necessarily prevents overlap in a 
highly-cluster region, such as Perseus, will struggle to find the true 
source population. Since one is normally interested in the material 
that will go on to form the final star, GAUSSCLUMPS would seem 
the more representative algorithm choice. 

Analysing the clumps found by two of the most popular 
clump-finding routines side-by-side, underlines the caution that one 
must exercise when comparing results from distinct studies on sim- 
ilar datasets in the same region. Different methods can produce very 
different results from the same data, providing clump populations 
that are not directly comparable. Future continuum and molecular- 
line surveys will need to carefully design their source extraction 
procedures to ensure objects are found at the desired length-scale 
without simply using an automated algorithm blindly. 
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APPENDIX A: CLUMP CATALOGUES 



This paper has been typeset from a Tj^X/ I^TbX file prepared by the 
author. 
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Table Al. Properties of SCUBA clumps found with CLFIND. The full version of this table is available as Supplementary Material to the online version of this 
article. 



Sub-region 


Clump 


RA-' 


Dec " 






5850 '' 






ID 


(J2000.0) 


(J2000.0) 


(pc) 


(Jy beam^ ' ) 


(Jy) 


(M.) 


NGC 1333 


1 


03:29:11.88 


+31 


13:09.4 


0.020 


3.47 


17.56 


21.02 


NGC 1333 


2 


03:29:03.23 


+31 


15:57.5 


0.030 


3.34 


22.23 


26.61 


NGC 1333 


3 


03:28:55.27 


+31 


14:36.5 


0.043 


2.76 


21.71 


25.98 


NGC 1333 


4 


03:29:10.01 


+31 


13:39.4 


0.024 


2.18 


10.71 


12.82 


NGC 1333 


5 


03:29:01.36 


+31 


20:30.5 


0.052 


1.40 


23.48 


28.11 



" Peak clump flux position. 

^ Clump effective radius, the geometric mean of the clump major and minor 'sizes' (\/Sizel x Size2). Each size is formed from the standard deviation of the 
clump co-ordinates about the centroid position, weighted by the pixel data values. 

Clump peak flux. 

Clump total flux. 

Clump mass derived from the total flux assuming optically thin, isothermal dust at 12 K (see section [33) . 

Table A2. Properties of SCUBA clumps found with GAUSSCLUMPS. The full version of this table is available as Supplementary Mateiial to the online version 
of this article. 



Sub-region 


Clump 


RA" 


Dec" 




r c 
'^peak 


^850 


^850 ' 




ID 


(J2000.0) 


(J2000.0) 


(pc) 


(Jy beam" ' ) 


(Jy) 


(M,5) 


NGC 1333 


1 


03:29:11.41 


+31:13:21.4 


0.021 


3.29 


18.87 


22.58 


NGC 1333 


2 


03:29:02.99 


+31:15:54.5 


0.017 


3.15 


12.25 


14.66 


NGC 1333 


3 


03:28:55.27 


+31:14:36.5 


0.013 


2.63 


6.82 


8.16 


NGC 1333 


4 


03:29:01.36 


+31:20:33.5 


0.021 


1.32 


9.40 


11.25 


NGC 1333 


5 


03:29:11.19 


+31:18:27.4 


0.016 


1.18 


4.87 


5.83 



" Peak clump flux position. 

Clump effective radius, the geometric mean of the clump major and minor 'sizes' (\/Sizel x Size2). Each size is formed from the standard deviation of the 
clump co-ordinates about the centroid position, weighted by the pixel data values. 

Clump peak flux. 
^ Clump total flux. 

Clump mass derived from the total flux assuming optically thin, isothermal dust at 12 K (see section [33) . 
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